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Abstract 

We extend the self-consistent Ornstein-Zernike approximation (SCOZA), first 
formulated in the context of liquid-state theory, to the study of the random 
field Ising model. Within the replica formalism, we treat the quenched ran- 
dom field as an annealed spin variable, thereby avoiding the usual average over 
the random field distribution. This allows to study the influence of the distri- 
bution on the phase diagram in finite dimensions. The thermodynamics and 
the correlation functions are obtained as solutions of a set a coupled partial 
differential equations with magnetization, temperature and disorder strength 
as independent variables. A preliminary analysis based on high-temperature 
and 1/d series expansions shows that the theory can predict accurately the 
dependence of the critical temperature on disorder strength for dimensions 
d > 4. For the bimodal distribution, we find a tricritical point which moves 
to weaker fields as the dimension is reduced. For the Gaussian distribution, 
a tricritical point may appear for d slightly above 4. 

Key words: Disordered systems, Ornstein-Zernike equations. Random Field 
Ising model. 
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Typeset using REVT^ 



I. INTRODUCTION 



Our understanding of the full phase diagram of randomly disordered magnetic systems 
in finite dimensions is severely hampered by the lack of a theory that takes into account 
the effect of fluctuations in an approximate but sensible fashion. This is probably the 
reason why the nature of the paramagnetic to ferromagnetic transition in the random field 
Ising model (RFIM) is still under debate after nearly twenty years of intensive studies 
(for reviews see [|T]-§1). The influence of the random field distribution on the order of the 
transition is one of the important open problems. Whereas mean-field theory which 
should be valid in high dimensions, predicts that some distributions give rise to a tricritical 
point, neither numerical simulations nor high-temperature series expansions [[10| -[T2 



have been able to yield a clear-cut answer yet. Recent numerical determinations of the 
ground states at zero temperature do not seem to clarify this issue |]13|,|T3[. In addition, 
as suggested by several analytical works [ITsHlSH, it may be that the true phase diagram is 



more complicated than anticipated, with the occurrence of an intermediate "glassy" phase 
signaled by a breaking of symmetry in the replica formalism. The fact that the RFIM, for 
a certain range of temperatures and fields, has a complicated energy landscape (so that the 
mean-field equations for the local magnetizations have many solutions pO[) also explains the 
failure of standard renormalization group perturbation theory pT| , p2| and the breakdown of 
dimensional reduction. 

In this work, we propose a theory for the RFIM that allows an approximate study of 
the influence of dimensionality on equilibrium properties. This is done by extending to this 
model the so-called self-consistent Ornstein-Zernike approximation (SCOZA) developed by 
Hoye and Stell for simple fluid and lattice-gas systems. Recently, the SCOZA has been 
shown to provide a very good description of the properties of the three-dimensional Ising 
model, even in the close vicinity of the critical point |^. We expect that the same approx- 
imation scheme will be useful in the case of disordered systems, as well. In a previous work 



25| , we studied as a first application the site-diluted Ising model, showing that this approach 
can indeed provide an accurate description of the dependence of the critical temperature 
on dilution. To overcome the lack of translational invariance of the Hamiltonian due to the 
random disorder, we applied the replica method in an unusual way, replacing the original 
system by n + 1 coupled systems with translationally invariant interactions. In the following, 
we shall use the same procedure and assume that there is no violation of replica symmetry. 
Our formalism can be generalized to study a possible replica symmetry breaking, but this 
question will be adressed in a forthcoming paper dealing also with the spin glass problem 



The paper is organized as follows. In section II, we derive the replica-symmetric Ornstein- 
Zernike equations for the RFIM. In section III, we introduce successively the Random Phase 
approximation (which is just another way of obtaining the mean- field thermodynamics) 
and the Optimized Random Phase approximation, which represents a first improvement on 
mean-field theory for hard-spin systems in finite dimensions. We then derive in section IV 
the SCOZA partial differential equations for the Gaussian and bimodal distributions of the 
random field. In section IV, we analyze the solution in terms of high-temperature and 1/d 
series expansions. A summary and a discussion are provided in section V. 
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II. REPLICA-SYMMETRIC ORNSTEIN-ZERNIKE EQUATIONS 



The RFIM is defined by the Hamiltonian 

n = -jY, - + h^)(yi (1) 



<i3> 



where J > 0, CTj = ±1, and < ij > indicates that the sum is over pairs of nearest neighbor 
sites of a d-dimensional lattice. H is a uniform magnetic field and the local fields hi are in- 
dependent random variables distributed according to some common probability distribution 
V{h). Quenched thermodynamic averages are defined by 



<A>T = Tr[exp{-(3n)A]/Tr[exp{-(3n)] (2) 

where (3 = l/{kBT) and the overbar denotes an average over the random field distribution. 
If not stated otherwise, we shall concentrate on the class of distributions that depend on a 
single positive parameter ho which measures the strength of disorder. They thus satisfy 

V{h)dh = V{T)dT (3) 

where r = h/ho and 'P(t) is independent of ho- This includes the much studied Gaussian 
distribution, 

Vih) = i2nhl)-'^' exp[-hy{2hl)] , (4) 
and the bimodal distribution, 



Vih) = ^[5{h~ho) + 5{h + ho)] . (5) 

In the first case, the quenched variables Tj are continuous (soft) spin variables whereas 
they only take the values ±1 (Ising-like) in the second case. The crucial point is that both 
quenched and annealed variables, Tj and at respectively, are present. Such quenched-annealed 
two-species systems have been previously introduced in the context of liquid-state theory 
to describe continuum fluids adsorbed in porous media [27,2S]. The replica method was 
used to obtain a set of exact equations for the pair correlation functions |^ and to derive 
thermodynamic relations The same technique can be used here. We first introduce n 



copies of the annealed spin variables and consider the Hamiltonian 

nn = -J J2 crtcr]-hoJ2^X-HY,^t ■ (6) 

<ij>,a i,a i,a 

Now, in contrast with what is usually done, we do not perform the average over the 
disorder variables to get an effective Hamiltonian (see however Appendix A). We treat 
the variables } and {rj} on an equal footing and we consider HrepiWi}, {^i}) = Ti-n — 
1/(3 InP(rj) as the Hamiltonian of a mixture of {n + 1) spin species. The corresponding 
free energy is J-'rep = —1//? InfTr exp(— /^Ti^ep)], where the trace is taken over {cr"} and {xj}. 
The average free energy of the RFIM is then given by 
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= lim - [exp{-(3J^rep) - 1] , (7) 

p n~^0 n 

and the average magnetization m = < dj >t is given by 

m = lim - < af >rev (8) 

a 

where < ... >rep denotes the average with respect to TLrep- The disconnected and 
connected correlation functions, Gdisij^) = < c"o >t< c"r >t — and Gconlr) = 
< (ToCTr >T — < ctq >r< >T, are related to the replica correlation functions C^^^r) =< 

CrSc^r >rep " < CTq >rep< CT^ >rep {a, b = l..n) hj 

= lim -r^-TT V ^"'(r) (9) 
n^o n(n - 1) ^ ^ ^ 



and 



Geo„(r) = hm - V G""(r) - G,,,(r) . (10) 



We also introduce the correlation functions Goo(i") = tot^—t"^ and 6*01(1") = tq < (Xr >t— 
which are related to the replica correlation functions (^^^(r) =< roTr >.,.ep — < tq >rep< 

Tr >,.ep and G^'^{r) =< TqC^ >rep - <ro >rep< K >rep by 



Goo(r) = limG""(r) (11) 



and 



Goi(r) = limi5^G°"(r). (12) 

a 

One has 6*00(1") = (t^ — r^) 5r,o for the distribution functions of uncorrelated random fields 
that are considered here (in principle, the case of correlated fields can be studied with the 
present formalism as well). Because of the hard spin condition cTj = ±1, one also has the 
sum-rule 

^^^(r = 0) = 1 -m^ (13) 

where (jii(r) = Gconij) + Gdisiy) (the Fourier transform of Gii(r) is the structure factor 
measured, e.g., in scattering experiments 0). Eq. (13) is equivalent to the so-called "core" 



condition in a lattice-gas On the other hand, 

G^onir = 0) = 1 - g (14) 



where q = < (Ji >^ is the standard spin-glass order parameter. 

We now provisionally assume that the external magnetic field is nonuniform, replica- 
dependent, and has an extra component that acts on the spin Tj. One has =< 
0",? >rep= —dTrepI dHI whcrc o uow represents not only the n replicas l...n but also species 
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and and m° stand for Tj and < Tj >rep, respectively. The Legendre transform that takes 
the fields i?f into defines the Gibbs potential 

n 

Qrep^:Frep+J2J2^>i (15) 

which satisfies 

- a< ■ ^^^^ 

Qrep generates the direct correlation functions (or proper vertices in field-theoretic language) 
which we define by 

Since C^j = — d^J-'rep/dHf^dHj, we have a set of Ornstein-Zernike equations 

n 

i c=0 

which, in the limit of a uniform replica-independent magnetic field, become in Fourier space 

n 

J2G^'{k)C'\k) = 6a,b . (19) 

c=0 

By taking the limit n —>■ and assuming replica symmetry, we finally obtain 

Goo(k) = (20a) 
(^00 (kj 



G'con(k) = - (20b) 

^cony'^) 

Gd^s{^) = [Si^ - C,U^)]^;z^, (20c) 

Coo(k) C'con(k) 

goi(k) = - . ^''^^^ , (20d) 

Coo(k)Ce.„(k) 

where the C's are related to the corresponding replica direct correlation functions by relations 
similar to Eqs. (9-12). Note that the first equation decouples from the other ones, as it should 
be, and that Gqi = Gio and Cqi = Cio by symmetry. These equations, hereafter called the 
replica-symmetric Ornstein-Zernike (RSOZ) equations, represent the starting point of our 
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study. Apart from a few notational changes, they are the same as those originally derived by 



Given and Stell [29| for a quenched-annealed mixture in the context of liquid-state theory (see 
also Refs. p^ , |33[] for an application to a lattice-gas model of a fluid in a disordered matrix). 
Because of the Legendre transform, m is now a control variable at our disposal instead of H, 
and in the following all quantities will be considered as functions of the three independent 
variables m, J = (3cJ (where c is the coordination number of the lattice) and = (3ho- We 
shall be especially concerned with the behavior of the susceptibility x(^) ^o) = dm/ d{(3H) 
given by 

X(m, J, ho) = Geon(k = 0) = ^ . (21) 

In the approximate theories that are discussed below, the divergence of x at fixed ho defines 
the spinodal line in the T — m plane. For a symmetric distribution (i.e., V{h) = V{—h)), 
and if the transition is continuous, the critical point is reached when the spinodal meets 
the magnetization curve at m = 0. Alternatively, we can locate the critical temperature by 
plotting x~^(m = 0) as a function of J at fixed ho ■ 



III. RANDOM PHASE AND OPTIMIZED RANDOM PHASE 

APPROXIMATIONS 

A. Random Phase Approximation 

When one turns off the exchange interaction in the RFIM Hamiltonian Ti, all quenched- 
averaged quantities can be calculated straightforwardly by direct averaging over the random 
field distribution. Hereafter, we shall call the system where J = the reference system. One 
then has 



m = tanh (3{hoT + H) (22) 

which can be inverted to get jSH as a function of m and ho (at least under the form of an 
infinite series). Then 



Qref = tanh^ p{hoT + H) (23) 
becomes a function of m and ho. For instance, 

qrefim, ho) = + (1 - m'^)'^[hl - 2(1 - 2m^)ho + \{81m^ - 90m^ + 17)hl + 0(/i^)] (24) 
for the Gaussian distribution, and 

qrefim, ho) = + (1 - m^)^[hl - ^(Sm^ + 1)^ + j-{225m^ + SOm' + 17)hl + 0{hl)] 

(25) 

for the bimodal distribution. 

In the reference system, all correlation functions are purely local. 
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GTlir) = (1 - qref) 5r,0 (26) 

G7/s{r) = {Qref-m')6.,o (27) 
Gll^{r) = Vref5r,o (28) 



where Vrefi'm, ho) = rtanh/3(/ioT + H) — rm. From the RSOZ equations, we then get 

C::ar) = 5r,o (29) 

-L — Qref 

C.7/(r) ^ ^.0 (30) 

I J- - Qref) 

Cll^{v) = -^r^L^ 5r,o (31) 

where Cqo = (t^ — t^)^^ (cqo = 1 for the symmetric Gaussian and bimodal distributions 
defined by Eqs. (4) and (5)). Note that a special feature of the bimodal distribution is that 
C^ifir) = 0, as can be readily checked. 

Once a reference system has been chosen, the Random Phase approximation (RPA) con- 
sists in adding to the direct correlation fonctions of the reference system the pair interactions 
which had been turned off (see, e.g., Ref. [0). In the present case, since there is no di- 
rect interaction between distinct replicas in the Hamiltonian Tirep (because the average over 
disorder has not been performed explicitly), the RPA in Fourier space writes 

C'r^(k)=C':i^(k)-JA(k) (32a) 
C'r^(k)=C',7/(k) (32b) 

C'or^(k) = (k) (32c) 

where A(k) = l/c^gexp(zk.e) is the characteristic function of the lattice and e denotes a 
vector from the origin to one of its nearest neighbors. Using Eqs. (21) and (29), this leads 
to 

1 Qref /no\ 

Xrpa = z (33) 

-1- — Zrpa 

where zrpa = </(l — Qref)- The same result is obtained by differentiating the mean-field 
expression of the magnetization 0] , 



m = tanh (3{cJm + /iqt + H) , (34) 
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with respect to /3H. Therefore, the above RPA description is equivalent to the standard 
mean- field theory of the RFIM [|,|]. Note however that the RPA free energy is identical to 
the mean-field free energy only when it is obtained by integration of the susceptibility. As 
will be stressed in section IV, there are several routes to obtain the thermodynamics from the 
correlation functions, and when the latter are only known approximately the different routes 
may not lead to the same results. For spin systems, the free energy or the Gibbs potential 
can be computed either from the susceptibility given by Eq. (21) or by integrating with 
respect to temperature the enthalpy which is itself expressed in terms of the pair correlation 
functions. For simplicity, we only discuss in this section the results of the susceptibility 
route. 

Then, for a symmetric distribution, and if the transition is second-order, the boundary 
between the ferromagnetic and paramagnetic phases is given in the RPA by zjiPA{fn = 
0,Jc,ho) = 1, i.e., 

Jc[l - qref{m = 0, ho)] = 1 . (35) 

In this approximation, the transition remains second-order along the whole phase boundary 
for the Gaussian distribution [Q. On the other hand, it becomes first-order at sufficiently 
large disorder strength when the random field distribution has a relative minimum at zero 
field [^. For instance, the tricritical point occurs at Jt = 3/2,tanh^(/io,i) = 1/3 for the 
bimodal distribution. Whether this scenario still holds in finite dimensions (i.e. when c is 
small) is an open question. 

B. Optimized Random Phase Approximation 

One shortcoming of the RPA is that the sum-rule, Eq. (13), is not satisfied. Indeed, 
from the RSOZ equations, we have 

G^f\k) = ^ ~ '^^'i + 'i^^f-'^^ (36) 

1 - ZRPA A(k) [1 - ZRPA m? 

(this expression leads to the well-known Lorentzian and Lorentzian-squared terms in the 



mean field structure factor near /c = g]). Introducing the lattice Green's function |35 

1 /-TT ik.r 

{2TiY J-n 1 - ^A(k) 



1 fTT ik.r 

P(r,z) = -— A (37) 



and going back to real space, we get 

Gfi^^(r) = (1 - qref)P{r, zrpa) + {qref - m2)[P(r, zrpa) + ZRPAP'{r, zrpa)] (38) 
where P'{r, z) = dP{r, z)/dz. We thus have 

Gfi^^(r = 0) = (1 - qref)P{zRPA) + {qref ' m^)[P{zRPA) + ZrpaP'{zrpa)] (39) 

where P{z) = P{r = 0,z). In general, this is different from 1 — m^. 
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To cure this problem, one may simply add to CHlijc) a state- and field- dependent per- 
turbation potential that is different from zero at r = and chosen in such a way that Eq. 
(13) is satisfied (of course, in the true system, the observables cannot depend on the value 
of the spin-spin interaction at r = 0). Using the terminology of hquid-state theory p4| , 
we call this approximation the Optimized Random Phase approximation (ORPA) ||3^. We 
thus write 

(^£r^(k) = Ce(m, J, h,) - J A(k) (40) 
whereas Eqs. (32b) and (32c) are left unchanged. Introducing zqrpa = J Icc, we now have 

GO^^^(k) = ^-l^^— + f^'f ~ (41) 

1 - ZORPA A(k) (1 - QrefY [1 - Zrpa A(k)]2 

and the sum-rule, Eq. (13), writes 

'-^P{zorpa) + 1;'''~'^' C-^?[P{zorpa) + zorpaP\zorpa)] = 1 - (42) 

which is viewed as an implicit equation for zorpa{^, J, ^o)- 

The susceptibility Xorpa = zorpa/IA'^ " zqrpa)] diverges when zqrpa = 1- For a 
symmetric distribution, and if the transition is second-order, the inverse critical temperature 
is then given by 

^(1) I qref{m = 0,ho) P{l) + P'{l) ^^3^ 
Jc [1 - qrefim = 0, ho)]^ 

In the absence of quenched disorder (i.e. when Hq = 0), one has g^e/ = and the inverse 
critical temperature is given by Jc = -P(l), i.e., Pc = PiX)/{.cJ). In this case, the ORPA 
is identical to the mean-spherical approximation (MSA) ||3^. (The MSA can be extended 
to the RFIM by choosing as reference the system where J = and = 0; then, Eqs. 
(32b) and (32c) are replaced by C^f^{k) = and C^^^{k) = -Jiq, which amounts to a 
linearization of the ORPA expressions with respect to hQ.) 

For d oo, P{z) ^ 1 and P'{z) 0, and we recover from Eq. (43) the mean- field 
equation for the critical temperature, Eq. (35). On the other hand, for d < A, P'{z) diverges 
at 2; = 1 (whereas P{1) diverges for d < 2), and the critical temperature, solution of Eq. 
(43), is driven to zero. This unfortunate drawback is shared by any Ornstein-Zernike theory 
that assumes that Ccon{T^) has the same range as the exchange interaction and that makes 
use of the susceptibility route. In this case, one finds from the RSOZ equations that the 
critical exponents rj and f], defined by Gcon{T^) ~ r^'^^^^'' and Gdis{T^) ~ r^'^^^^*' when r — 00 
at the critical point, are both zero. As a consequence, there is no critical point at nonzero 
temperature for (i < 4 because this would lead to the unacceptable result that G'ii(r) does 
not decrease to zero at long distances. 

By using the expansion of P{z) about z = 1 for d > 2 and d 7^ 4, 6, 8... [^ , 



Piz) = [P(l) + 61(1 -z) + 62(1 - zf + ...] + (1 - zY^^-'lco + Ci(l -z)+ C2(l - z)^ + ...] 

(44) 
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(for d = 4,6,8, there are logarithmic corrections), it is easily found from Eq. (42) that 
for 4 < d < 6 the ORPA gives rise to the critical exponents of the random-field spherical 
model I^H], i. e. z/ = l/((i — 4), 7 = 2/((i — 4) and 6 = d/{d — 4). The mean-field exponents 



are recovered for d > 6. Therefore the ORPA susceptibility route leads to a d ^ d — 2 
dimensional reduction that is a direct consequence of the Ornstein-Zernike approximation 
for the direct correlation functions. On the other hand, it is easy to show that the enthalpy 
route yields classical (mean-field) critical exponents. Both routes give back the mean-field 
results in the limit d ^ 00 [a more complete study of the ORPA for hard-spin systems with 
quenched disorder will be presented elsewhere |3D|). 



An illustration of the predictions of the ORPA is given in Tables I and II for the 5- 
d hypercubic lattice (c = 10, P(l) = 1.156308, P'(l) = 0.778633). The inverse critical 
temperatures obtained from Eq. (43) are compared to the best available estimates obtained 
from a fifteen-term high-temperature series expansion of the susceptibility [0 (the disorder 
strength is measured here in terms of the reduced variable g = h'l/{cJY = h^/J"^; this 
differs from the g defined in Ref. [|1^ by the factor 1/c^). Although the ORPA represents a 



clear improvement on mean-field theory, we see that the agreement with the "exact" results 
deteriorates significantly as the field increases. For the bimodal distribution, Eq. (43) has 
no solution when g > 0.1199 and one may suspect that the transition becomes first-order 
for the highest values of the field. In order to decide on the existence of a tricritical point, 
we expand the inverse susceptibility around m = along the critical isotherm and look for 
the change of sign of / I 1"^"^ at m = and J = Jc- Since Xorpa ~ (1 — zorpa) when 
zoRPA 1~, this amounts to expanding the solution of Eq. (42). Using the expansion 
given in Eq. (44), we find after some calculations that the two conditions for tricriticality 
are Eq. (43) and 

^[P(l)+P'(l)] = -1 (45) 

where v42(/io) is the coefficient of in the expansion of {qref — m^)/(l — QrefY around 
m = 0, 



AM-^^^^^^^^^^^^ ^ (46) 

^ ^ (l-t2)5 (1-^2)2 ^ ^ 

where t = tanh(r/;,o). For the Gaussian distribution, Eqs. (43) and (45) have no solution and 
no tricritical point appears in finite dimension. On the other hand, there is always a solution 
for the bimodal distribution, and one finds for instance g^^^^ = 0.117 and J^^^^ = 1.85 
for the 5-d hypercubic lattice. By comparing to the mean- field predictions, g^^^ = 0.193 
and Jf = 1.5, we see that there is a range of fields where the transition is driven first 
order by fluctuations, in agreement with the conclusion of a previous high-temperature series 
analysis [|lT|. For more general symmetric random-field distributions, a careful analysis of 
Eqs. (43) and (45) shows that the conclusion of Aharony based on mean-field theory 
remains unchanged for the ORPA: a sufficient condition for the occurence of a tricritical 
point is that V{h) has a minimum at zero field. 
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IV. SELF-CONSISTENT ORNSTEIN-ZERNIKE APPROXIMATION 



As is well known in liquid-state theory and has been mentionned above, solving the 
Ornstein-Zernike equations with an approximate expression of the direct correlation func- 
tions like in the RPA, the ORPA or any other approximate closure relation such as the 



Percus-Yevick approximation or the hypernetted chain equation ||3J], generally leads to 
thermodynamic inconsistencies. In the language of magnetic systems, this means that dif- 
ferent Gibbs potentials are obtained depending on whether one uses the susceptibility or 
the enthalpy routes (the enthalpy is defined by £^ = ^ -f- TS = U + MH where S and lA 
are the entropy and internal energy, respectively, and M = Nm is the total magnetiza- 
tion). For instance, in the pure Ising model, the former route corresponds to the double 
integration of the equation = C(k = 0) = d'^ {f3Q / N) / dm? with respect to m (at 
constant T), and the later to the integration of the Gibbs-Duhem relation £ = djSQ/djS 
with respect to f3 (at constant m). The requirement that the two routes lead to the same 
results is thus embodied in the relation dC{\i. = 0)/d(3 = d"^ {S / N) / drr? , provided that 
the appropriate initial conditions are satisfied. Since the enthalpy per spin is given by 
S/N = —J/NJ2<ij> < '^i^j >T= —cJ/2[G{y = e) +m^], this relation, together with the 
Ornstein-Zernike equation, may be viewed as an exact equation for the pair correlation 
function. 

The extension to the RFIM is straightforward. Let us consider the general variation of 
the average free energy due to variations of the three control parameters J, and H. We 
have 

dJ^ = —dJ < aiaj >t — dho < (Xj >t — dH < cTj >t ■ (47) 

<ij> i i 

Then, from the definitions of Gii(r) and G'oi(r), we readily find that 

d{l3G/N) = -^(Gii(r = e) + m^)dJ - (Goi(r = 0) + fm)dho + (3Hdm . (48) 
On the other hand, we have from Eq. (17) 

4.(k.0,^«, (49) 



SO that we get the three "Maxwell relations" : 



a4on(k = 0) ^ ia2Gn(r = e) 
Yj = "'"2 d^- ' 



gG'oi(r = 0) ^ 1 .9G'n(r = e) 

dJ 2 dhn ^ 



dConjk = 0) ^ g^Goi(r = 0) 
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Since the enthalpy density is given from the Hamiltonian 7i by 



8/N = -y [Gn(r = e) + m^] - h,[G,,{v = 0) + fm] , (53) 



the combination of Eqs. (50) and (52) yields 

80,^^ = 0) d\£/N) 



(54) 



dp dm? 

which generalizes the self- consistency relation for the pure Ising model discussed above and 



used in Ref. 



The SCOZA strategy is now the following. We assume that the three direct correlation 
functions Cconiy)^ Cdisir) and Coi(r) have the same range as in the RPA but with values (for 
the present problem at r = and/or r = e) which are state- and field-dependent. We thus 
write in Fourier space 

^^^^^^^(k) = ce(m, J, ho) [1 - z{m, J, ho)\m (55a) 
^^.^^^^(k) = c,(m, J, ho) (55b) 

C'oT''^^(k)=coi(m,J,/io) (55c) 

(for notational simplicity, we do not use the index SCOZA for z,Cc,Cd and Cqi). It follows 
from the RSOZ equations that 

G!^n'''^ir) = -Pir,z) (56a) 

Cc 

G'l^'^ir) = '^^^^^§-JzP{r, z)] (56b) 
G^^cozA^r) = -^P{r,z). (56c) 

Cc Coo 

We then impose that the exact relations, Eq. (13) and Eqs. (50-52), be satisfied. This leads 
to a set of partial differential equations (PDE) in the unknown quantities z,Cc,Cd and cqi. 
This is not enough, however, because only two of Eqs. (50-52) are independent provided 
that the appropriate initial conditions are satisfied (for instance Eqs. (50) and (51)). We 
thus need an additional relationship between the pair correlation functions. 

This extra equation is readily found in the case of the Gaussian probability distribution 
which has the special property that / A TV{T)dT = f{dA/dT)V{T)dT. Choosing A =< 
Gi >T yields < cxj >t tj = d < ai >t /drj = [3ho < aiaj >t — < ai >t< o-j >t, and one 
gets the exact relationship 

Goi(r) = ho ao„(r) . (57) 
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Since cqo = 1, this gives, when inserted in the RSOZ equations, 

C'oi(k) = -h^ . (58) 

We thus take cqi = —ho in the SCOZA equation (55c) and we use the sum-rule, Eq. (13), 
to express q as a function of Cc, 

~2 P{z) - Ce(l - m^) 

= ^° + P(.)+.n.) • 

Using the fact that the Green's function at nearest-neighbor separation satisfies P(r = 
e, z) = [P{z) — V\lz, we finally obtain from Eqs. (50) and (51) two coupled PDE's in 
z{m, J, ho) and Cc{m, J, ho), 

Ar^i = _^ri ^(^)-^ + (1 - - ^) 1 (60b) 

cj 9(/.g)^c, z cJp(z) + zP'(z)J- ^'''^ 

When setting = P{z)/{1 — m?) in these two equations, the first one reduces to the PDE 
for the pure nearest-neighbor lattice-gas {ho = 0, with the usual substitution p = (1 + m)/2 



and w = A J) which has been studied in Ref. |2^; see also the discussion in section V. In 
fact, the exact relationship, Eq. (57), suggests that it was unnecessary to introduce the 
correlation function G'oi(r) in the Gaussian case and that one could have averaged over 
the quenched disorder from the outset, as is usually done. We show in Appendix A that 
the same expressions of the correlation functions and of the thermodynamic quantities are 
obtained within the SCOZA when one uses this alternative route. 

The case of the bimodal distribution is somewhat more complicated. As already noticed, 
it also has a special property, namely, C^if{r) = 0. However, Cdisir) is not zero when J 7^ 
and there is no reason a priori to set Q(m, J, ho) = in Eq. (55b). The solution to this 
problem consists in introducing an additional independent variable that allows to control 
the mean value of the random-field. Indeed, Eq. (17) tells us that 



^ij' (^Om'^dm'i ^^^^ 

^ J 



which gives, in the hmit n — » and for a uniform magnetic field, 



In consequence, we derive from Eq. (48) two additional Maxwell relations 

dCoi{\i = {)) _ l d^Gii{T = e) 
aJ ~ 2 df dm 

9(701(1^ = 0) ^ 92Goi(r = 0) 



(63) 



dho ^ df dm ' ^^^^ 
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Only three of Eqs. (50-52) and Eqs. (63-64) are independent but we now have with Eq. (13) 
the right number of equations to determine unambiguously 2;, Cc, Q, and Cqi. In Appendix 
B, we apply this procedure to the asymmetric bimodal probability distribution V{h) = 
p6{h — ho) + {1 — p)6{h + ho) . Since f — 2p—l, the parameter p can be varied independently 
of m, J and ho to change f. This leads to three PDE's, Eqs. (B3). Of course, it is 
significantly more difficult to solve three PDE's than only two and it is highly desirable to 
simplify the problem, especially if one is only interested in the case p = 1/2. Fortunately, it 
turns out that q = is a very good approximation to the full solution of Eqs. (B3) (more 
precisely, we show in Appendix B that = through order (3^ in the high-temperature 
series expansion of the solution of Eqs. (B3)). This is probably related to the fact that 
one has C^j^/(r) = even when p ^ 1/2. With this approximation, we are left with only 
three unknown quantities, z,Cc and cqi, which can be determined by using Eq. (13) and 
Eqs. (50-51). Returning to the case p = 1/2 and introducing f{m,J,ho) = — Cqi/cc as a 
new unknown variable, we eliminate Cc from Eq. (13) (cf. Eq. (B2) with r — f), 



Piz) 



m? - P[P{z) + zP'{z) 



(65) 



and we finally obtain from Eqs. (50) and (51) two coupled PDE's in z{m,J,ho) and 
/(m, J, ho) (cf. Eqs. (B3a) and (B3b) with r = /), 



dj'l-m^-p[P(z)+zP'{z)Y 2dw?^' J \ \ J' V y/j ^p^^^ 

+ fP'{z)} (66a) 

^ifPiz)] = -m'- f{P{z) + zP'iz))]^^^ + fP'iz)} (66b) 
dJ 2 dho zF{z) 



V. SERIES EXPANSIONS OF THE SCOZA SOLUTIONS FOR THE GAUSSIAN 

AND BIMODAL DISTRIBUTIONS 

To integrate the coupled PDE's, Eqs. (60) or (66), the initial values of 2;, or / for J = 
must be known. The quantity z is related to the second moment correlation length ^ defined 
by qe(k) = Gc(k = 0)[1 - i^k'^ + ...]. Specifically, one has z = c/{c + T^)- Therefore z = 
ioT J — since all correlation functions are local in the reference system (^ = 0). On the 
other hand, Cc,re/ = 1/(1 ~ Qref) and in the bimodal case fref is given by = qref — tn^- 
As discussed in section III. A, g^e/ can be obtained by direct averaging over the random- 
field distribution. It is much more convenient, however, to calculate the properties of the 
reference system from the third "Maxwell relation", Eq. (52). When J = 0, this equation 
leads to 

d Ccrefjm, ho) r d'^ 1 /^^x 

dho Cc,ref{m,ho) 

for the Gaussian distribution, and 
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_^ 1 ^ d'^fre/jm, hp) 

for the bimodal distribution. These equations are solved by imposing the boudary conditions 
at m = : qref{0,ho) = l/y^27r) J^^exp{—u'^ /2) tanh^{uho)du for the Gaussian distribu- 
tion and gre/(0, /io) = tanh^(/io) for the bimodal distribution (morever, /^e/ and Cc,re/ are 
even functions of m). This is how we have obtained the series for qref in Eqs. (24) and (25). 

It may be noticed that Eqs. (60) and (66) have a different behavior in the limit 0. 
Indeed, z and Cc are even functions of ho whereas / is an odd function. Therefore, the solution 
of Eq. (66b) when /iq — > is / = 0, which leads to Cc = P{z)/{1 — rn?) so that z is identical 
to the solution of the equation for the pure system that has been considered in Ref. ||24|| . 
On the contrary, in the Gaussian case, the solution of Eq. (60b) is not Cc = P{z)/{1 — m^) 
when /lo 0, and one does not recover the pure system results. The discrepancy, however, 
is extremely small. 

In Ref. [Q, the numerical integration of the SCOZA partial differential equation for 
the pure Ising model was performed by rewritting this equation as a quasi-linear diffusion 
equation for which implicit predictor-corrector algorithms are available in the literature. 
The numerical integration of the coupled PDE's, Eqs. (60) or (66), is more difficult and we 
defer this task to a later work. Note however that an interesting feature of the theory is 
that a single run of integration steps sweeps the whole parameter space, so that the phase 
diagram in the T — ho plane can be obtained at once. In what follows, we only give some 
preliminary results obtained from series expansions. 

The SCOZA equations are indeed quite suitable for deriving high-temperature series. 
Since z vanishes when J — > 0, one can expand the Green's function P{z) in powers of z, 

P{z) = l + J2Pkz'' (69) 

k>2 

(with P2 = 1/c), and substitute into the PDE's. We then express z{m, J, ho),Cc{m, J, ho) 
and /(m, J, ho) as triple series in J, m and ho, 

z{m, J, ho) = (1 - m^) Zijkhfm^^J" (70) 



cAm 



J, ho) 



E 



ijk'''0 



(71) 



/(m, J, ho) = J2 kkhf^'rn'^ J' , (72) 

by using the fact that z vanishes when m = ±1. Eventually, we must gather in Cconi)^ = 0) 
all terms at a given order in (3 in order to obtain from Eqs. (21) and (55a) the high- 
temperature series expansion of the zero-field susceptibility. A careful analysis of the PDE's 
shows that to calculate x(^ = 0) through the order /?" it is sufficient to consider the triple 
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series, Eqs. (70)-(72), up to finite values of k : imax = n, jmax = n — i and kmax = n — i 
or {n — i)/2 for the Gaussian or bimodal distribution, respectively (this is for Eqs. (70) 
and (71); in Eq. (72), kmax = n — i — 1/2 oi (n — i — l)/2). The crux of the calculation 
is that the coupled PDE's reduce at each order in J, m, to a system of hnear algebraic 
equations in the unknown coefficients Zijk and Cc,ijk or Zijk and fijk- In consequence, the 
whole calculation can be performed with reasonable effort using a symbolic computation 
software like MAPLE or MATHEMATICA. The results for the hypercubic lattice in general 
dimension are 

Xim = 0) = 1 + 2rf(/3J) + {-2d + AS - g){^Jf + (4rf/3 - + - Mg){l3Jf 
+ (lOrf/3 + 16c/V3 - 24c/3 + IQdt" + (4c/ - 12^^)^ + 2g^){pjf 

+ (-28C//3 + 88c/V5 + 2Ad? - 64rf^ + 32c/^ + (-8rf/3 + 24^^ - ?,2d^)g + lMg'^){(3Jf + ... 

(73) 

for the Gaussian distribution, and 

Xim = 0) = 1 + 2d{pJ) + {-2d + 4rf2 _ g){l3jf + (4rf/3 - 8^^ + - Ug){l3Jf 
+ (10c//3 + 16c/V3 - 24rf^ + 16c/^ + (4rf - 12^^)^ + 2//3)(/?J)^ 

+ (-28(^/3 + 88rfV5 + 2Ad^ - 64rf^ + 32rf^ + (-8c//3 + 24^^ - ?,2d^)g + 14rfc/V3)(/5^) 

(74) 



for the bimodal one (here g = Ad'^g as in Ref. |1T2|]). Comparison with the exact series |T2 
shows that the SCOZA series are exact through the fourth-order term (at the order jS^, the 
only inexact coefficients are those of d and d"^ which take the values -28/3 and 88/5 instead 
of -116/15 and 16, respectively). More generally, the numerical values of the higher-order 
coefficients are in remarkable agreement with the exact ones. As an illustration, the terms 
of order f3^^ are given in Appendix C for both distributions (note that at each order in /3 
the coefficient of the highest-order term in g is exact: this is because the reference system 
is treated exactly as a boundary condition to the PDE's). Similar series expansions can be 
obtained for the structure factor G'ii(k = 0) and for other types of lattice (one only has to 
use the corresponding Green's function in the PDE's). 

We now turn to more quantitative predictions. Assuming that = 0) ~ (1 — 
T/Tc)'"' near T = T^, we have performed a Dlog Pade analysis of the high-T series (since 
we hope to have a numerical solution of the PDE's in the near future, we have not tried to 
apply more sophisticated methods of analysis such as the ones used in Ref. [0). For the 
critical temperature, the convergence between the different approximants is fairly good and 
the results in c? = 5 using a [5/6] approximant are given in Tables I and II (this approximant 
has been chosen because of its smooth behavior as a function of g). We see that the SCOZA 
predictions are in very good agreement with the estimates obtained from the exact high- 



temperature series expansion |]12[. The phase diagrams displayed in Figures 1 and 2 show 
that the improvement on the RPA and the ORPA is indeed remarkable. On the other hand, 
the predictions for 7 are rather sensitive to the order of the approximant and are therefore 
less reliable than for the critical temperature. 7 must anyhow be interpreted here as an 
effective exponent since we expect to obtain the random-field spherical-model exponents 
asymptotically (this is not yet proved, however). The increase of 7 with g has been noted in 
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preceding studies [0 but for a serious discussion we must wait until the numerical solution 
of the PDE's is available. 

To discuss the possible occurence of a tricritical point in hypercubic lattices, it is more 
convenient to consider 1/d expansions. These can be obtained either by using the l/d expan- 
sion of the Green's function, P{z) = 1 + {2d) + 3z^/ (2(i)^ + ..., or by using the conventional 
scaling J = 1/ {2d) to reorganize the J-expansions. In the pure Ising model, this latter pro- 
cedure readily yields the 1/d expansion of the critical temperature around the mean- field 
d ^ oo limit. In the RFIM, the calculation is more complicated because the mean-field 
critical temperature is itself solution of an implicit equation, Eq. (35). In consequence, we 
cannot use the conventional high-temperature /3-expansion as a starting point, but rather 
the J-expansion at constant ho = [3ho. In order to get the systematic corrections around 

mean- field theory, it is more convenient to take qref{fn = 0) = tanh^(^o''") instead of Iiq as 
independent variable in the PDE's. All this complicates the formal procedure, especially 
in the Gaussian case, and we have computed the expansion only through order 2. For the 
bimodal distribution, we find that the inverse zero-field susceptibility is given by: 

1. s 1 ~ J"^ (l + x)(l-2x) .P I ,2, . 

-J{l%x^ - 32x3 ^ ^2x^ + a; _ 1)] + (75) 

where x = qrefi'm = 0, ho) = tanh^(/io). The coefficient of the 1/d term can be shown to be 
exact. The above expression yields the expansion of Jc around the mean- field inverse critical 
temperature J^^, solution of Eq. (35). The result is reproduced here only to first order: 

jMF ^ 2d 1-x^^P 1-x^'P -2{gx^F)y^ ^d'' ^ ' 

where x^^ = 1 — 1/J^^^. The tricritical point is obtained by solving simultaneously the 
equations x~^(m = 0) = and d'^x~^ /^m? |m=o= (since one always has classical mean-field 
critical exponents in a l/rf expansion, we take here 5 = 3). This yields 

and 

~ o.^^^ 0.1320 0.4950 ^,1, 

,, = 0.1927-^^ ^ + P8) 

(recall that g = h^/J"^). We thus approximately locate the tricritical point in c/ = 5 at 
Jt = 1.771 and gt = 0.147 (see Fig. 2). This confirms the conclusion reached with the 
ORPA: there is a range of field strengths (for instance, 0.147 < g < 0.193 in 5-d) where 
the transition is driven first-order by fiuctuations. The anomalous behavior of 7 observed 
at large g in Table II is probably related to a crossover to this tricritical behavior. 
For the Gaussian distribution , we find 
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, I z J^ xK(x) , .P ,2^^,,, K(x)-2il-xf 

X-\m = Q) = - J- [i-x + —^] + ——{-K{x){l+x) ^ ^ 
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l-x 2d' (l-x)2J (2rf)2'3 ' ' 

+^l(r^*-^(--) + + + (1 - + 0(1) (79) 

where x = 1/^/121^) exp(— ^2/2) tanh^(-u/io)(iM and K{x) = {2ho)^^ dx/dfiQ. The func- 
tion ii'(x) can be calculated once for all; one has -ft'(x) = l—Ax + 9x'^~2Ax^ + 96x'^ + 0{x^), 
and since x = 1 - v/(2/7r)/io^ + 0(/io^) when /iq ^ oo 0, ~ (7r/4)(l - x)^ when 

x ^ 1. Here too, the expansion is exact to first order in 1/d. The resulting expansion of Jc 
is: 

1.^11. ,/'^''-'::s.j .o(i) (80) 



JMF 2d' {l-x^^Pf -2gK{x^^Py 'rf2' 

The expansion given by Eq. (79) also permits to calculate the expansion of g* , the 
variance of the random-field distribution for which Tc = 0. We find 

^•-!-^<^ + ^)-(2^4-^^)^0'^'- 

Since there is no tricritical point in infinite dimension, we cannot perfom a perturbative 
expansion of Jt or gt around mean-field theory. However, we can expand the equations 
X~^(m = 0) = and d'^x~'^ldvn? |m=o= in powers of l/rf and look for a solution at a 
given value of d. For > 5 there is no indication that a tricritical point occurs. On the 
other hand, the calculation at order l/rf^ suggests that the transition becomes first order 
for sufficiently strong disorder when d < 4.2. This is certainly a rough estimation which 
must be confirmed by the numerical solution of the PDE's. The fact that d ^ 4 might be 
a "critical" dimension for the Gaussian RFIM has been also suggested by preceding high-T 
series analysis [|rT[ . 



Finally, we have also calculated the high-T and 1/d series expansion of the quantity A 
defined as 

A = li,„ 1 G^Ok^ ^ l^C|A^-0) _ c^^^ ^ 



It has been argued in the literature ^2] , ^ that A should be equal to 1. In agreement with 



the analysis of Ref. ||T2|, we find that A is always finite (this is here a direct consequence 
of the OZ approximation for the correlation functions) and close, but not exactly equal to 
unity. 



VI. CONCLUSION 

The calculations of the preceding section, based on the analysis of high-T and 1/d series 
expansions, shows that the thermodynamically self-consistent Ornstein-Zernike approxima- 
tion (SCOZA) accurately predicts the critical temperature Tc{hQ) of the RFIM in dimension 
d > 4 for both the Gaussian and bimodal distributions (at least in the range of disorder 
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strength where we can compare to the available "exact" results). For the bimodal distri- 
bution, we find that the phase transition becomes first-order for sufficiently strong random 
fields, in agreement with mean-field theory 0, and that the tricritical point moves to weaker 
fields as the dimension is reduced. A tricritical point may also appear in the Gaussian case 
for dimensions slightly above 4, but the numerical solution of the partial differential equa- 
tions is needed to settle this question. This numerical solution, when available, will also 
permit to probe the critical region and to calculate the effective exponents above the critical 
point (the analytical determination of the true asymptotic exponents, especially on the co- 
existence curve, and the elucidation of the scaling behavior of the theory are also challenging 
tasks). The accuracy of the theory shows that the assumption that the connected and dis- 
connected direct correlation functions have the same range as in the RPA (i.e., essentially 
the same range as the interaction potentials) is quite reasonable, as long as one does not 
approach the critical point too closely or go to very low temperatures (the exact structure of 
the correlation functions is certainly more complicated at low temperatures, as illustrated by 
the exact results in 1-d |^^). The main shortcoming of this Ornstein-Zernike assumption is 
that 7] = f] = which forces Tc to be zero for d < A. On the other hand, the SCOZA can be 
easily generalized to the n-component version of the RFIM, as it has been done for the pure 
Ising model [^. One can then show that the theory becomes exact in the spherical-model 
limit n —>■ oo. This means that the SCOZA should be even more accurate for the classical 
random-field X-Y or Heisenberg models. Finally, it must be stressed that all the above 
results assume that replica symmetry is not broken in finite dimensions. This assumption 
must be justified (or invalidated) by studying the stability of the replica symmetric solution. 
This will be done in a forthcoming paper devoted also to the application of the SCOZA to 
the Ed wards- Anderson spin-glass model [p^ . 
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APPENDIX A: 



In the case of the Gaussian distribution, one usually averages exp(— /37Y„) over disorder 
to get an effective Hamiltonian 

<ij>,a i,a,b i,a 

The species does not appear any more but the replicas are now interacting. The corre- 
sponding RSOZ equations are 

Gco„(k) = (A2a) 



Working first in replica space and then taking the limit n — > 0, it is easily shown that 
the differential expression for the Gibbs potential is now 

d{pg/N) = -^(Gii(r = e) + m^)dJ - hoGcon{r = 0)dho + PHdm (A3) 

(f = for the Gaussian distribution). Because of Eq. (57), one thus gets the same 
Maxwell relations as Eqs. (50-52). On the other hand, some care must be taken in 
computing the internal energy or the enthalpy. The internal energy is obtained from 
U. = lim^^o dUeff /On where Ueff satisfies the Gibbs-Duhem relation Ueff = d{(3J-'eff)/d(3. 
Since the Hamiltonian Heff is temperature-dependent, one has the unusual expression 
Ueff —< Ti-eff + P dHeff/dp >eff where < .. >e// denotes the average with respect to 
Heff- Prom Eq. (Al), this gives 

^ = - Y E[^""(r = ^) + - (^^0 E[^"'(r = 0) + ^'^^"i - ^ E ^" ■ 

a a,b a 

Taking the limit n — > and using Eqs. (10) and (12), one obtains 

S/N = -y [G'n(r = e) + m^] - phlG^onir = 0) , (A5) 

which, owing to Eq. (57), reduces to Eq. (53), as it should be. 

Since the interaction between replicas is on-site, the range of the SCOZA direct correla- 
tion functions Cf,^*^-^^(r) and C|j^*^-^^(r) is also unchanged, 

C^T^''{k) = I ho)[l - z'{m, J, ho)X{k)] (A6a) 
<^fr^^(k) = c:,(m,J,M, (A6b) 
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so that 



Gern^^^(r) = >(r,.') (A7a) 

Gf^^^^(r) = A[.'P(r, .')] . (A7b) 

It is then clear that the requirement that Eqs. (13), (50) and (51), together with the exact 
initial conditions at J = 0, be satisfied implies that z' — z and = Cc and d^ — 0^ — c^i- 

APPENDIX B: 

We consider the asymmetric bimodal distribution 

V{h)^p5{h-ho) + {l-p)6{h + ho) . (Bl) 

One has f = 2p — 1 and Cqo = [4p(l — p)]"^. Introducing f{m,J,ho) = —Cqi/Cc and 
r^(m, J, ho) = (cqi/coo — Cd)/cl as new unknown variables, we use Eq. (13) to eliminate Cc, 
i.e., 

" l-m? -r''[P{z) + zP'{z)] ^ ' 

and wc obtain from Eqs. (50), (51), (56), and (63) three coupled PDE's in z{'m, J,ho), 
f{m, J, ho) and r(m, J, ho), 

^[ l-..^-r^[P(.)+.P-(.)] ] - - '-2d^^^^' - - - ^^^^^^ + ^^))]44r 

+ r''P'{z)} (B3a) 

4p(l - p)i[/P(^)] = ^^{[1 - - r^(P(^) + zP'{z))]^^^ + r2p'(^)} (B3b) 

(B3c) 

As explained in section V, the solution of the SCOZA partial differential equations can 
be easily obtained in terms of high-temperature series expansions. The remarkable feature 
of the solution of Eqs. (B3) is that one has q/c^ = 4p(l — p)p — = through order 
More precisely, for the hypercubic lattice in general dimension, we find 

Cd/cl = p2(l - pffm{l - m^f[-hm + Ug^'^{p - ^)]/?^ + 0(/3^) (B4) 

5 2 
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Note that this term cancels out when m = (this is no longer true for the next-order 
term). This implies that the zero- field susceptibility for the symmetric bimodal distribution 
is the same through order (3^ when calculated from Eqs. (B3) or from Eqs. (66). A more 
complete study based on a twelve-term series expansion shows that q = is a very good 
approximation when p = 1/2. On the other hand, it would also be worth studying also the 
full solution of Eqs. (B3). Indeed, it is easily seen that the same equations are obtained 
for the probability distribution V{h) = p6{h — ho) + (1 — p)6{h + hi) (just define r by 
h = l/2[(/io + hi)T + {Jiq — /ii]). In mean-field theory, this asymmetric RFIM has a rich 
phase diagram with several coexistence surfaces PT| and it has been claimed to be relevant 
to the liquid- vapor coexistence boundary of ^He confined in silica aerogel. 



APPENDIX C: 



To illustrate the accuracy of the SCOZA, we show here the terms a [15] in the high-T 
series expansion of the the zero-field susceptibility, x(^ = 0) = Yl^=o'^[^]{l^'^Y ■ These 
results can be compared to the exact ones given in Ref. |T^. One finds 

a[15]= mmd}^ - 229376^1'* + 1835008(i^V3 - 3923968rf^V5 + 19566592^^745 

+62568448rf^7315 + 2393632768ciV5775 + 1309641268736ciV225225 + 344496180736^78775 
-250795760924032(^7259875 + 1073182266768952832^7212837625 
-2633146526329112672^7212837625 + 687607165382389072c/742567525 
-6687777748364066392(i7638512875 + 522635516549738848rf/212837625 
+ (-114688d^^ + 638976^^2 - 1306624rf" + 3549184d^73 - 6071296ci7l5 - 26430464^745 
-5725256704^710395 + 530937363712c/7289575 + 87605199818176^72027025 
-24080214624448^7135135 + 1584841452347248^76081075 - 75721982136512c/7405405 
+357823952018272d/6081075)^+ (184320rf^^ - 768000ci^° + 1136640rf^ - 3053568^75 
+2969344^715 - 4845952^7105 - 447235513408^751975 + 2060640064288^751975 
-11034060991936^7155925 + 1826181780848(i731185 - 4442599992^/275)^^ 
+ (-547840^73 + 509952ci^ - 4678144^79 + 57472^^ - 10716992^7315 
+3003435424^7945 - 8105781856(^7945 + 21474530792^72835 - 4915534784^/1701)^^ 
+ (384256^73 - 187776ci^ + 1368512ci79 + 2739616^745 - 85830176c/7l35 
+262363504c^7315 - 10055056ci/105)^^ + (-355008^75 + 52160^73 - 666224^715 
+23348344c/7225 - 14219344c^/225)/ + (1560464^745 + 1000808^745 - 199312(i/9)c/6 
-4709644cic/7315 

(CI) 

for the Gaussian distribution function, and 

a[15]= 32768rf^^ - 229376^^^ + 1835008rf^73 - 3923968d^75 + 3940352d^79 + 59213824^^7315 
+7336961024ci7l7325 + 261515720192^745045 + 26451169582592c/7675675 
-1086932402689792^71126125 + 27530806058894912^75457375 
-877599448749913984^770945875 + 255326320300912^715795 
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-2202231489060446384^7212837625 + 1549572126103322224d/638512875 
+(-114688^^=^ + 638976ci^2 - 1306624^^^ + 3528704d^73 - 1995776^75 
-8554496(i7l5 - 167141888^7297+14174567936^76435 + 10235242804672^7225225 

-29834386736192^7155925 + 596984424197072^72027025 
-407009452618592^72027025 + 310791876465968^/6081075)^ 
+ (151552rf" - 1853440^^73 + 892928^^ - 1246720^73 + 29252864(i7315 
-3510912^735 - 65936770624ci77425 + 1238778408016^731185 
-11162124651232^7155925 + 7782983079632^7155925 - 19615564589086^/467775)^^ 
+ (-846848^79 + 3586048^715 - 6394880^727 - 1122176^745 - 1242688^7105 
+997169632^7405 - 52631828384^78505 + 4259397904^7675 
-817609382692^/127575)^^ + (1751296^763 - 433024^715 + 1367232^735 
+115455808^74725 - 1154097376^74725 + 9855809872^733075 - 20089162036^/99225)/ 
+ (-16538048^74725 - 474176^7567 - 85810544^7 14175 + 543224408^770875 
-3984248rf/6075)^^ + (65701616^7467775 + 18472712(i7l55925 + 23549048^/1403325)^^ 
-35397196d^742567525 

(C2) 

for the bimodal distribution. 
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l3cJ 


7 


9 


RPA 


ORPA 


SCOZA 


"exact" 


SCOZA 


"exact" 





0.1000 


0.1156 


0.1139 


0.1139 


1.031 


1 


0.08 


0.1087 


0.1370 


0.1311 


0.1315 


1.094 


1.12 


0.10 


0.1112 


0.1438 


0.1369 


0.1369 


1.142 


1.13 


0.12 


0.1138 


0.1514 


0.1425 


0.1428 


1.133 


1.14 


0.14 


0.1165 


0.1601 


0.1492 


0.1493 


1.148 


1.142 


0.15 


0.1179 


0.1649 


0.1528 


0.1529 


1.158 


1.144 


0.18 


0.1224 


0.1818 


0.1656 


0.1665 


1.192 


1.245 


0.20 


0.1257 


0.1959 


0.1759 


0.177 


1.219 


1.28 


0.25 


0.1349 


0.2489 


0.2090 


0.213 


1.190 


1.30 



Table 1: Inverse critical temperature and critical exponent 7 for the 5-d liypercubic 
lattice (Gaussian distribution). The SCOZA predictions have been obtained from a [5/6] 
Dlog Pade analysis of the high-T series expansion of the zero- field susceptibility. The "exact" 
results are taken from Ref. [Q, except for g = which is taken from Ref. 
in Ref. corresponds to our g multiplied by = 100. 



Note that g 
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7 


9 


RPA 


ORPA 


SCOZA 


"exact" 


SCOZA 


"exact" 





0.1000 


0.1156 


0.1139 


0.1139 


1.031 


1 


0.05 


0.1057 


0.1296 


0.1251 


0.1253 


1.077 


1.0 


0.07 


0.1085 


0.1379 


0.1312 


0.1314 


1.087 


1.1 


0.08 


0.1100 


0.1430 


0.1347 


0.1349 


1.090 


1.11 


0.09 


0.1116 


0.1494 


0.1386 


0.1389 


1.086 


1.11 


0.10 


0.1134 


0.1577 


0.1431 


0.1435 


1.077 


1.135 


0.11 


0.1154 


0.1700 


0.1473 


0.1488 


1.005 


1.15 


0.12 


0.1175 




0.1546 


0.1552 


1.073 


1.18 



Table 2: Same as Table 1 but for the bimodal distribution. 
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FIG. 1. Phase diagram of the random field Ising model with a Gaussian distribution in d = 5. 
Comparison of the mean-field (long-dashed line), ORPA (short-dashed line) and SCOZA (solid 
line) predictions to the "exact" phase boundary (dots). The SCOZA results are obtained from 
the high-temperature series expansion of the zero-field susceptibility. In principle, the SCOZA 
line should continue down to = but no reliable Fade approximant has been found at low 
temperatures. 
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FIG. 2. Same as Fig. 1, but for the bimodal distribution. Only the portion of the curves 
corresponding to a second-order transition is represented. The crosses locate the tricritical points. 
The SCOZA tricritical point is estimated from a l/d series expansion at second order. 
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